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^ '. Abstract 

The gravitational field of matter and that of antimatter could differ. This might be one signature 
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of quantum gravity. We show that primordial Big Bang Nucleosynthesis restricts such a possibility. 
PACS numbers: 



The symmetry between matter and antimatter arises naturally in conventional quantum 
field theory. With the priors of Lorentz invariance, hermiticity, and locality, the theory is 
CPT conserving. Among other consequences, CPT conservation implies that the inertial 
masses of a particle and its corresponding antiparticle are equal. When technically feasible, 
this identity among masses has been checked in the laboratory. General relativity is based 
on the equality of inertial and gravitational masses, and this hypothesis should apply in 
principle both to matter and to antimatter. However, it is not at all obvious that this 
assumption will be tenable when going from classical general relativity to a quantum theory 
of gravitation. The path to quantum gravity is not yet established, but one thing we know 
for sure is that we have to go beyond conventional field theory; it follows then that we cannot 
invoke the kind of matter- antimatter symmetry we alluded to. In fact, there is the claim 
ll LJ lii L| about the possibility of a fundamental asymmetry between matter and antimatter 
that could manifest as a signature of quantum gravity (for a review see |5[ and references 
therein) . 

It is not a trivial issue to probe the gravitational properties of antimatter and compare 
them with those of matter. Apart from our own work, that we will present below, we are 

n n 

aware of only one instance where this comparison has been possibly done p, l2|- in the 
supernova 1987A collapse, a huge flux of neutrinos and antineutrinos was emitted. The flux, 
when crossing the Earth, left a signal in underground detectors. The statistical analysis 
concluded that ja|, at the 90% confidence level or more, at least one v e e event was detected, 
the rest being u e p. Assuming v e and v e were indeed detected, and taking into account 
that all events were recorded in about a 10 second period, one can put a restriction on the 
geodesic deviation of these particles when following a path from the supernova in the Large 
Magellanic Cloud (LMC) until the Earth. This is expressed conveniently in terms of the 
post-Newtonian parameter 7 (7 = 1 for all particles and antiparticles in Einstein's theory). 
One gets |fi b| 

It* -7pJ<KT 5 -10- 6 (1) 

The bound (JTJ) on the (weak) equivalence principle for matter and antimatter is subject to 
some caveats. As we said, it is necessary that at least one of the detected events corresponds 
to a v e . In addition, if the violation of the equivalence principle is due to new long-range 
forces, then the limit (0) is relaxed if the range is shorter than the LMC-Solar System 
distance. Also, screening effects might invalidate (JTJ), as is discussed in |9(. 



We should also say that, after the successful antihydrogen production |lCJ,|ri||, a labora- 
tory limit on conceivable differences between the fall of hydrogen and antihydrogen in the 
Earth gravitational field will be in principle possible [12J. 

This last experimental proposal and the limit <$J$ are concerned with the effects of a 
gravitational field on matter and antimatter. We would also like to compare the gravita- 
tional field produced by bulk matter with the corresponding gravitational field produced 
by antimatter. This is very hard of course because we need a large macroscopic sample 
of antimatter and measure the gravitational field it produces. Fortunately, the very early 
universe offers a "laboratory" for such a test. To start with, one has comparable densities of 
matter and antimatter. Also, the particle and antiparticle energy densities are responsible 
for the gravitational Hubble expansion. In the standard Big Bang model, the expansion rate 
of the universe, measured by the Hubble parameter H, is given by the Friedmann equation 
in terms of the energy densities p i: 

H 2 = yGnJ2p* ( 2 ) 

i 

where G^ is the Newton's gravitational constant. In (J2J), the contribution of curvature and 

cosmological constant has been neglected. This is justified, since we shall apply the equation 

in the very early universe, specifically in the period of Big Bang Nucleosynthesis (BBN). 

We see that in (J2J) the densities pi of the different species i present in the early universe 

contribute to the expansion with the same strength Gn, so that, as a consequence, matter 

and antimatter enter with the same gravitational constant. The root for that is found in 

Einstein equations 

W v _ 1 gl *, R = g7r q nT w (3) 

The Friedmann equation (J2J) can be deduced from Einstein equations, if we use the 
Robert son- Walker metric to calculate the Ricci tensor R^ u and the Ricci scalar R, and 
describe the energy-momentum tensor T^ u as a perfect fluid. The energy densities pi are 
contained as a perfect fluid source in T^ v and they contribute with the same strength G^. 
If matter and antimatter produce different gravitational fields, we would expect a viola- 
tion of Einstein equations ©, and as a consequence the Friedmann equation (J2J) would also 
be invalid. Even if we do not have a fundamental theory in this case, we shall proceed in 
the following way to constrain such matter/antimatter asymmetry. 



We will make use of the fact that the Friedmann equation gives a fair description |l3j| of 
the universe expansion in the BBN period. It is reasonable to think that if the Friedmann 
equation (J2J) is to be altered, the modifications are bound to be not excessively drastic. 

With this argument in mind, we shall postulate a minimally modified Friedmann equation 
where matter and antimatter enter with different strengths, 

H 2 = y [G N p M + G N {l + 8 1 )p 1 + G n (1 + 6 a )pa] (4) 

Several comments are in order. We have left Gn multiplying the energy density matter con- 
tribution pm, since G^ is the usual magnitude measured for the gravitational field produced 
by a source of matter. The energy density of electrons and neutrinos will contribute to pu 
(protons and neutrons are present in the BBN era but, in practise, it is the relativistic par- 
ticles that dominate pu)- A modification appears for antimatter, which enters our equation 
(HJ) with a strength Ga = Gn(1 + 5a)- We will have positrons and antineutrinos contributing 
to the antimatter energy density pa- Also, there are photons with an energy density p 7 that 
contribute to the universe expansion. Once we admit the possibility of breaking the mat- 
ter/antimatter symmetry, we should let the photons enter the modified Friedmann equation 
with yet another gravitational constant G 7 = Gat(1 + 5 7 ) different in principle from both 
matter and antimatter. 

We shall apply BBN arguments to constrain l A and K BBN has already quite a long 
list of limits to different models |TJ] . Let us first summarise how it works in our case. 

In the standard BBN, it is crucial when weak interactions freeze out, since from this 
moment on the interconversion between protons and neutrons is ineffective, and primordial 
nucleosynthesis starts. A modification of the expansion rate H makes the freeze out temper- 
ature different from the standard one and this affects the prediction of the primordial yields 
of light elements. We will calculate the theoretical yields using the modified Friedmann 
equation Q with the parameters 6a and 5 7 . Then we will compare the theoretical predic- 
tions with observational data and require agreement. In this way we shall obtain bounds on 
the parameters 5a and 5 1 . For the predictions we need to specify the number of neutrinos 
N v and the neutrino/antineutrino asymmetry which is described by the chemical potential 
\x v . This is easy to understand: changing the number of neutrinos changes the expansion 
of the universe, on the one hand, and, on the other hand, p v enters in the proton/neutron 
interconversion reaction rates. We shall assume the values of the standard BBN: N„ = 3, as 



collider measurements tell us, and \x v = 0. In the standard BBN analysis one gets that [x v /T 
has to be small 15], with the electron neutrino asymmetry more restricted that the others. 
Neutrino oscillations seem to imply that the limits for the electron neutrino apply also to 
the other neutrinos 16]. Even if we work in non-standard BBN, we think it is reasonable to 
set \x v = in our analysis, and in this way we do not have an extra parameter to play with. 
In addition, the predictions depend on the baryon number density. So, still another 
parameter 77 which is this density normalised to the photon number density is needed in 
our calculation. Until very recently BBN was the only way we had to get information on 
1). This situation has changed with the observational results (l2j on the cosmic microwave 
background, leading to the allowed range 

r} = (6.13 ±0.25) x 1(T 10 (5) 

We will use (|HJ) as an input. Finally, the primordial yields of light elements also depend on 
the neutron lifetime; we take r n = 885.7 ± 0.8 sec J18j |. 

The success of BBN in constraining models is of course not only due to the theoretical 
framework we have but also to the fact that we possess observational data. The data is not 
without debate since some hypotheses have to be done to go from the observed abundances 
to the primordial abundances. We adopt the recent measurements coming from |l9(], where 
the deuterium detection in quasar absorption systems at high redshift leads to the following 
ratio of deuterium to hydrogen by number 

D/H = 2.78+^xl0~ 5 (6) 

Unfortunately, the determinations of primordial 4 He suffer from discrepancies between 
the different estimates. We prefer to take a conservative attitude and, instead of choos- 
ing a specific measurement among all, we adopt a range which encompasses the existing 
measurements. The observational primordial 4 He mass fraction is taken to be in the range 

0.228 <Y P < 0.248 (7) 

However we will comment below what would happen if we follow a more recent analysis. 

We are now ready to explain how we actually get our results. We run a modified version 
of the numerical Kawano code 20[. We first introduce as inputs the central values of i] 
and r n that we mentioned before. For each value of 5a and 5 7 we get a prediction for the 



abundances of deuterium and helium-4. Comparing the predictions with the data in © 
and (J7J) we get two allowed bands. These are shown in Fig. QJ Requiring consistency for 
deuterium leads to the band between the two solid lines, while the dashed lines enclose the 
allowed band for 4 He. It is no surprise to get a band since there is some cancellation between 
the effects of 5 a and 5 1 if the sign of these two parameters turns out to be opposite. The 
bands have not the same slope and thus we get an allowed region when we consider both 
elements, leading to the bounds —1.4 < 5a < 0.2, —0.2 < 5 7 < 1.2. 

We have arrived at our main objective: to show that BBN provides a limit on the possi- 
bility that the gravitational field of antimatter and radiation differs from the gravitational 
field of matter. We got limitations on 5a and 5 1 for rj = 6.13 x 10~ 10 . We should now allow 77 
to vary in the experimental error range and strictly speaking we should do that introducing 
7] as a statistical variable. We will not do it this way because, as we said, the observational 
4 He, which is one of the main inputs, is taken with a large error to take into account mea- 
surements that are not consistent at one a among themselves. We hope the situation will 
improve in the future, but for the time being we do not undertake the task of treating r\ 
statistically. For our purposes, we limit ourselves to find the allowed region in the 5a, 5^ 
plane for rj = 5.88 x 1CT 10 and r] = 6.38 x 1CT 10 , namely the one-a observational extremes 
coming from WM AP. We keep r n fixed since this does not introduce any appreciable change. 
The results are shown in Figs. |2(a)| and |2(b)| We conservatively take as the allowed range 
for our parameters the union of the ranges in the figures, which is 

- 1-6 < 5 A < 0.4 

-0.4 < 5 1 < 1.5 (8) 

This is the main numerical result of our work. 

In our analysis, we have used a 4 He abundance that is large enough to comprise different 
measurements. In fact, there are recent analysis leading to values Yp that are not consistent 



with both the measured D/H © and rj (JSJ). For example, taking Yp = 0.238 ± 0.005 21 1 
leads to this situation. Allowing a non-zero 5a and 5 1 would make a consistent scenario 
where Yp, D/H and r] agree. We think however that the tension between the observables is 
likely to be solved by an improvement of the observations. So we regard our work as giving 
a constraint on antimatter gravity, rather than giving a solution to a crisis. 

We finally would like to stress that if in our modification of the Friedmann equation Q 



we have some relation between 5a and <5 7 then the bounds will be tighter. For example, 
assume the following modified Friedmann equation 



r 2 ™ 



H 2 = — [G(l + 5)p M + G Pl + G(l - 5) PA ] (9) 

Here matter and antimatter couple differently but the departure from the radiation coupling 
is somehow symmetric. This is what would be expected in theories with a "graviphoton" 
22 1 and has been shown to be consistent [23] with energy- conservation arguments 24]. In 
(JHJ) we have to identify G^ = G(l + 5). The limits from BBN on 5 are 

-0.02 < <5 < 0.12 (10) 

The bound is now tighter because now we have © that has one free parameter instead of 
(0J) with two parameters. We would have another one-parameter model if we put <5 7 = and 
let only vary 5 a- The bound would then be 

-0.52 <5 A < 0.05 (<5 7 = 0) (11) 

The reverse situation is that 5a = 0, 5 y ^ 0. Physically, it amounts to say that antimatter 
and matter behave identically while it is the gravitational field of radiation that is different. 
Notice that an alternative physical interpretation is that fermions and bosons produce a 
different gravitational field. In fact, this possibility has been considered in [25], where the 
authors use BBN to constrain this fermion/boson asymmetric model. 
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FIG. 1: Constraints on the (5a, S~,) plane coming from the requirement that the predicted abun- 
dances are consistent with the observed values. We show the constraints from deuterium as a solid 
line, and from 4 He as a dashed line. We have put r\ = 6.13 x 1CP 10 . 
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FIG. 2: Region in the (<5,4,<5 7 ) plane where the modified Friedmann equation leads to abundances 
in agreement with observation. We display the limits coming from deuterium (solid line), and 
4 He (dashed line). The values for the baryon number density are (a) rj = 5.88 x 10~ 10 , and (b) 
r} = 6.38 x 1(T 10 . 
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